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Abstract 


The  release  of  diatomic  species,  in  particular  AlO,  into  the  upper  atmosphere 
at  twilight  results  in  luminous  clouds  that  display  the  resonance  electronic - 
vibrational -rotational  spectrum  of  the  released  species.  The  intensity  distribution 
of  the  vibrational-rotational  bands  is  temperature  dependent.  Assuming  thermal 
equilibrium  of  the  released  species  with  ambient,  the  ambient  temperature  can  be 
determined.  Using  the  AlO  molecular  parameters,  Franck-Condon  factors,  and 
the  solar  intensity  in  the  upper  atmosphere,  the  population  of  AlO  in  the  vibrational 
and  rotational  levels  and  the  relative  band  intensities  are  calculated,  tabulated,  and 
charted.  These  charts  will  be  useful  for  upper  atmosphere  temperature  measure¬ 
ments  from  AlO  luminous  clouds. 


Foreword 


This  report  was  prepared  for  the  Air  Force  Cambridge  Research  Laboratories 
by  Ove  Harang.  The  work  was  initiated  under  Contract  No.  AF  1 9 (628) -5 7 31  with 
Northeastern  University,  Boston,  Massachusetts,  and  completed  under  Contract 
No.  AF  61  (052)-906  with  the  Auroral  Observatory,  Tromso,  Norway. 

Because  of  its  importance  as  a  guideline  for  upper  atmosphere  temperature 
measurements  from  fluorescent  chemical  releases,  it  is  being  published  in  an 
AFCRL  series  to  give  it  wider  distribution. 
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AlO  Resonant  Spectrum  for  Upper  Atmosphere 

Temperature  Determination 


I.  IM  ROHM  THIN 

In  recent  years,  rocket  release.-,  of  TMA  (trimethyl  aluminum)  trails  were 
performed  for  upper  atmosphere  wind  studies.  When  the  TMA  is  released  into  the 
atmosphere,  it  reacts  with  atomic  oxygen  in  a  chemiluminescent  reaction  forming 
AlO  (aluminum -monoxide).  If  the  trail  is  exposed  to  sunlight,  a  brilliant  sped  rum 
is  obtained  of  the  blue-green  system  of  AlO  (A^S*  -X**— +)  (Blanont  et  al,  li»t>l; 
Armstrong,  ltu*J;  and  Rosenberg  ct  al,  1  !>**•€). 

Another  way  of  creating  AlO  is  to  release  an  exploding  charge  seeded  with 
metallic  aluminum  into  the  uoper  atmosphere  (-Wilder  et  al,  I9t»,2,  19(d;  Arm¬ 
strong,  It**  I;  and  Ha  rang,  I’-*M). 

The  vibrational -rotational  band  intensity  distribution  of  the  AlO  spectrum  is 
temperature  dependent.  The  temperature  of  the  upper  atmosphere  can  be  deter¬ 
mined  through  analysis  of  the  AlO  spectrum  (Authier  ct  al,  lftoi,  IMJ,  I9M; 

Ha  rang.  it**  4;  and  .in! hie r,  l!**>4). 

1.  fill:  \N>  ROM  1 1 1.1 

Figure  1  shows  the  term  diagram  of  the  AlO  molecule,  with  transitions  indicated. 


(Received  for  publication  H  April  l !'"-*-) 


Figure  1 .  Term  Diagram  of  the  AlO 
Molecule 


The  only  electronic  transition  recorded  in  rocket  releases  is  the  "blue-green" 
system  A^2+  -X^£+.  The  ultraviolet  system  B2«-A2£  is  in  a  spectral  region 
where  photoelectric  sensitivity  drops  appreciably.  The  intensity  of  the  ultraviolet 
system  is  expected  to  be  weak,  because  the  solar  radiation  intensity  at  3000  A  is 

o 

one  fifth  that  at  5000  A:  if  the  f-values  of  the  two  transitions  are  the  same,  the 
absorption  intensity  of  the  blue-green  system  should  be  five  times  stronger  than 
the  ultraviolet  system. 

The  red  system,  which  shows  up  only  together  with  the  ultraviolet,  is  expected 
to  be  far  weaker  than  the  ultraviolet  because  of  the  v3  factor  in  the  Einstein  coeffi¬ 
cient.  In  the  present  model,  only  the  blue-green  system  is  considered  in  resonance 
radiation,  neglecting  possible  filling  up  of  the  A22+  level  from  the  level. 

Figure  l  shows  the  two  electronic -vibrational  levels  and  the  vibrational- 
potential  curves  plotted  as  a  function  of  internuclear  distance.  Each  vibrational 
level  is  built  up  of  rotational  levels.  The  position  of  a  certain  elect  ronic- 
vihrational -rotational  level,  T,  is  expressed  by  the  quantum  numbers  of  the 
electronic,  vibrational,  and  rotational  state. 


Figure  2.  The  A2S+  and  X2S+  Levels 

T  (n,  v,  K)  =  Te  (n)  +  Ty  (n,  v)  +  TK(n,v,K) 
where 

Ty(n,v)  =  wn(v  +  */2)  +  wxn(v+ 1/2)2  +  uyn(v+l/2)2 
TK(n,v,K)  =  B^KiK+l)  +  DnK2(K2+l) 

Bv  =  Bn  -  *(v  +  1/2) 

D„  =  4B*/Wn2. 

The  molecular  parameters  Te,  u,  uy,  B0,  and  a  are  taken  from  a  survey 
of  the  AlO  molecule  bv  Tyte  and  Nicholls  (lf64). 
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In  this  report,  the  following  set  of  parameters  was  used: 


State 

T 

A  e 

u 

ux 

Wy 

Be 

a 

A2S+ 

20.68895 

870.05 

3.52 

0.0 

0.6040 

0.0044  cm-* 

X2S+ 

0.0 

979.23 

6.97 

0.0 

0.6413 

0.0058  cm-1 

3.  THE  POPULATION  OF  v"  LEVELS  IN  THE  GROUND  STATE 

Following  the  formation  of  AlO  molecules,  they  will  collide  with  air  molecules, 
and  the  populations  of  the  different  rotational  and  vibrational  levels  will  be  given  by 
the  Boltzmann  distribution  law. 

For  vibration  levels 

Nv„  =  Qy'1  exp  -  (Gq  (vm  )/0.  6958T) , 
where 

Qv  =  2  exp  -  (Go(v")/0.6958T) , 


and 


G0  (v")  =  (u,  -  wx)  v»  -  «x(v")2. 
For  rotational  levels 


NKn  =  (2  K"  +  1)Qk-1  exp  -  (F(K")/0.6958T, 
where 


Qk 


2  (2  K"  +  1)  exp  -  <F(K")/0.6958T, 

K" 


and 


FK  =  B  K"  (K"4  1)  -  D(K")2  (K"  +  l)2  . 


Figures  3  and  4  show  plots  of  the  populations  Nv«  and  N^n  as  a  function  of 
the  levels  and  for  different  temperatures.  Here  Nvn  and  are  normalized, 


S  N 
v»  v" 


S 

K" 


1. 


1.0 


’,,t.  ■•••' _ r'-i-.r’/i  /i5j _ j^l 

1000  2000  3000 

e»*' 


Figure  3.  Population  of  Vibrational  Levels  of  AlO 

4.  SOLAR  INTENSITY 

To  compute  the  absorption  rates,  we  need  to  know  the  photon  intensity  of  the 
solar  radiation  on  the  top  of  the  atmosphere.  Figure  5  shows  the  intensity  meas¬ 
ured  in  photons  cm"^  sec'l/(cm"*)  and  in  photons  cm'^  sec**  Hz"*  as  a  function 

of  wave  numbers.  The  values  are  deduced  from  the  Handbook  of  Geophysics  (1960) 

0  0 

integrated  over  intervals  of  50  A  and  100  A,  including  the  Fraunhofer  lines  that  are 

o 

about  0.1  to  0.  5  A  wide. 

The  rotational  lines  of  AlO  are  about  0. 5  A  apart,  except  at  the  bandhead. 

For  vibrational  analysis,  the  entire  band  intensity  is  measured,  covering 

o  « 

about  1 0  to  30  A.  Solar  intensities  averaged  over  20  to  30  A  do  not  de%liate  much 

o 

from  values  averaged  over  50  A. 


Figure  4.  Rotational  Population  of  the  Ground 
Vibrational  Level  v"  =  0  of  AlO 


The  emission  from  a  certain  v'  -level  is  caused  by  absorption  from  several 
v" -levels,  so  this  will  tend  to  smooth  out  the  deviation  from  the  averaged  solar 
intensity  at  the  wavelength  of  an  individual  transition. 

There  exists  one  strong  and  broad  Fraunhofer  line,  H^,  which  happens  to  lay 
just  on  the  P-branch  of  the  0-0  transition.  For  600°K  this  correction  amounts  to 
less  than  10  percent  of  the  branch,  which  is  less  than  5  percent  for  the  entire  0-0 
transition. 

For  rotational  analysis,  it  is  of  importance  to  consider  the  Fraunhofer  corrcc 
tion  of  solar  intensity  lor  those  rotational  lines  which  are  strong  in  absorption. 
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Figure  5.  Photon  Intensity  on  the  Top  of  the  Atmosphere 


5.  THE  ABSORPTION-EMISSION  PROCESS 

In  order  to  compute  the  theoretical  intensities  of  bands  and  lines  of  an  assem¬ 
bly  of  AlO  molecules  in  the  upper  atmosphere  exposed  to  solar  radiation,  the  fol¬ 
lowing  model  is  considered, 

A  short  time  after  formation  of  AlO,  the  molecules  will  attain  thermodynamic 
equilibrium  with  the  atmosphere,  and  the  vibrational  and  rotational  levels  will  have 
a  population  according  to  the  Boltzmann  distribution  law. 

The  molecules  undergo  a  resonant  absorption  followed  by  emission.  The  tran¬ 
sition  probability  is  in  the  order  of  10^  -108  sec"*. 

The  emission -absorption  processes  will  tend  to  perturb  the  thermal  distribu¬ 
tion  of  the  lower  levels.  It  is  assumed  that  the  rate  of  the  absorption -emission 
process  is  relatively  slow,  so  that  there  is  a  sufficient  time  lapse  between  each 


absorption -emission  process  per  molecule  to  allow  the  molecules  to  readjust  to 
the  Boltzmann  distribution  through  collisions  with  atmospheric  molecules. 

This  depends  on  the  f -value  of  the  electronic  transition  of  the  AlO  molecule, 
the  collision  frequency  in  the  upper  atmosphere,  and  the  cross  section  for  colli- 
sional  deactivation  of  vibrational  levels. 

The  rate  of  excitation  from  level  X  to  level  A  is  expressed  by: 

=  NX  JLid  fxA  ev  =  NX  3-s±  ev  v 

dt  "me  v  3  he  gx 

where  9  is  the  solar  photon  intensity,  which  is  about  400  photons  cm*2  sec**  Hz"1 
in  the  region  where  AlO  absorbs.  The  f -value  is  unknown,  but  it  must  be  reason¬ 
ably  large  because  the  trail  is  almost  as  brilliant  as  a  sodium  trail.  The  density 
of  sodium  atoms  in  such  trails  is  probably  of  the  same  order  of  magnitude  as  that 
of  AlO,  and  the  f -value  of  sodium  is  0.98. 

If  fjr^  is  in  the  order  of  0. 1  to  1,  there  will  be  about  1  to  0. 1  absorptions  per 
second  per  molecule,  which  leaves  the  molecule  about  1  to  10  sec  to  attain  thermal 
distribution  between  each  excitation. 

The  number  of  collisions  required  for  thermal  redistribution  varies  within 
large  limits.  The  collision  frequency  in  the  upper  atmosphere  is  2  to  3  x  103  sec* 
at  100  km,  20  collisions  sec"1  at  150  km,  and  10  collisions  sec*1  at  170  km,  ac¬ 
cording  to  the  U.S.  Standard  Atmosphere  1962.  This  problem  is  further  discussed 
under  Section  1 1 . 

The  line  strength  S(A,X)  for  an  allowed  transition  is  defined  as: 

S  (A,X)  =  =  |<  Alplx  >\* 

where  p  here  is  the  dipole  moment.  The  degeneracy  gx  of  the  lower  state  of  AlO 
equals  two  (doublet). 

The  electronic  levels  are  composed  of  vibrational  levels,  and  the  expression 
for  absorption  from  a  lower  level  v'  to  an  upper  level  v"  is 

dNv,  a  .3 

=  Nx  N„  -5-—  w(V,  v")  S(V,  V)  0<v\  v»). 
d  t  a  v"  3  he 

Here  the  band  strength  is  (Tyte  and  Nicholls,  1964): 

S(v\  v")  x  Re2  0 -0.46  r(v',  v"))2  q(v',  v"). 

The  vibrational  levels  are  in  turn  composed  of  rotational  levels.  The  rate  of 
population  of  an  upper  rotational  level  K'  from  the  lower  level  K"  is 


— 
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dN 


d,:  *  NxNv  nk"  *'■  *"•  K")  S(V,  V)  ¥&■/;>. 

The  transition  considered  for  AlO  is  A2L+-X22+,  thus  K'  -  K"  can  only  have 
values  of  +  1  or  -  1 . 


K' -  K"  =  -1  (P-branch),  S(K\  K")  =  K'  +1  =  K" 

K’-  K"  =  +1  (R-branch),  S(K',  K")  =  K*  =  K"  + 1 

In  emission,  the  rate  of  electronic  transitions  from  A  to  X  is: 


dt 


=  na  aax  *  na 


.4  3 

^1*  -fts  S(A*X) 


3  he* 


«A 


Here  the  degeneracy  g^  also  equals  two. 
The  rate  of  vibrational  transition  is: 


=  N»  N  ,  -^.”4  v")3  s(V.  v")  , 

dt  A  3hc3 


and  a  similar  expression  for  a  rotation  transition  is: 


dNv* 


£  =  NA  N,„  Nk.  64"4v(V.  K'  V,  K-}3  s(v.>v>)  SiKVKl 
t  *  3hc3  2K’  +  I 


6.  THE  VIBRATIONAL  SPECTRIM 


The  population  rate  for  level  v'  caused  by  absorption  is: 


dNv, 
dt  V 


Nv«  S(v',  v")  v(v',  v")  0{v\  v«). 

V" 


The  depopulation  rate  caused  by  emission  from  an  upper  level  v'  is: 

dNV  ccNv.  £  v(v\  v")3  S(v\  v"). 
dt  ” 


These  two  rates  are  equal  in  tynamic  equilibrium  and  we  can  deduce  the  relative 
population  as 
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Nv,  oc 


5,  N  „  S(v',  v")  v  (v\  v»')  0(v\  v") 


*'<v\  v")3  S(v',  v") 


Figure  6  shows  the  plots  of  the  relative  population  Nv,/Nvi  _  q  as  a  function  of 
temperature. 

Having  computed  the  values  of  Nv,/Nv,  _  q  ,  the  intensities  of  the  bands  can  be 
computed: 


.  hi  v"  Nv"  S*v'*  V^v'*  v")0(v'.  v")  , 

*',v”)ccX. - y~_  - — - - S(V,  v")  v(v',  v")3  . 

S(v',  v")  v(v',  v")3 


Figure  6.  Population  of  the  Upper  States  for  Various  Temperatures 
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Figure  7  gives  the  relative  intensities  of  bands  within  the  Av  =  -1,0,  and  +  1 
sequence  as  a  function  of  temperature.  The  sequence  which  is  the  most  sensitive 
to  temperature  variation  is  the  +  1  sequence.  Unfortunately,  the  intensity  of  this 
sequence  is  the  weakest  of  the  three  sequences  (see  Figure  8).  The  strongest  one, 
the  Av  =  0  sequence,  shows  almost  no  temperature  dependence. 


I 

8 


>• 


£ 


Figure  7.  Intensity  Ratio  of  Rands  Within  Different 
Sequences  of  AlO  as  a  Function  of  Temperature 


Figure  8.  Intensity  Ratio  of  Sequences  as  a  Function 
of  Temperature 


i\ tknsitiks ok  tutu.  >kqi  km*:* 


So  fir,  the  band  intensities  within  each  sequence  and  their  temperature  varia¬ 
tion  were  deduced.  There  is  also  a  temperature  dependence  of  the  total  intensity 
ratio  of  different  sequences. 

I'sing  the  intensities  computed  in  Section  S  ail  the  intensities  within  each 
sequence  can  be  added.  The  intensity  ratios  Av  *  -t/Av  *  0  and  Av  *  I/Av  *  0 
are  plotted  in  Figure  8  as  a  function  of  temperature. 

The  advantage  of  measuring  the  sequence  intensities  is  that  is  it  not  necessary 
to  use  a  high -resolution  spectrometer.  A  photometer  with  interference  filters  of 
SO  A  half -width  will  suffice. 
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The  optical  throughput  of  a  photometer  can  be  made  very  big  compared  to  a 
spectrometer,  and  the  intensities  can  be  measured  with  high  precision. 

In  order  to  deduce  temperatures,  the  filter  profile  must  be  known  as  well  as 
the  transmission  of  the  optical  system  and  the  sensitivity  of  the  detector.  But  even 
without  precise  knowledge  of  these  constants,  we  can  record  variations  along  trails 
and  obtain  a  temperature  profile,  although  an  absolute  value  of  the  temperature 
might  be  difficult  to  obtain. 

It  is  noticed  that  there  is  very  little  temperature  dependence  for  temperatures 
below  500* K. 


I.  ERRORS  INVOLVED  IN  THE  C0HPLT4TI0N 

The  accuracy  of  the  computed  intensities  depends  mainly  on  the  following  items: 

1.  The  validity  of  the  model  adopted,  that  is,  complete  thermalization  of  the 
lower  v"  -levels  between  each  excitation -emission  event,  and 

2.  The  accuracy  of  the  parameters  used  in  the  computation. 

The  parameters  used  are  as  follows.  The  transition  strength  S(v',  v")  ex¬ 
pressed  by  the  product  of  the  Franck -Condon  factor  q(v*,  v")  and  the  electronic 
transition  factor  (1-0. 46  (v*,  v*))^,  and  the  intensity  of  the  solar  radiation  cor¬ 
rected  for  the  Fraunhofer  absorption  lines  0(v*,  v*). 

Until  recently,  all  temperature  deductions  were  based  on  the  approximation 
that  the  transition  strength  S(v',  v*)  is  expressed  by  the  Franck -Condon  factor 
q(v\  v")  without  the  r -centroid  correction. 

The  five  curves  in  Figure  9  illustrate  how  the  relative  intensities  of  bands  are 
affected  by  some  adjustment  of  the  transition  strengths  and  solar  intensities. 

Curve  I.  The  transition  strengths  S(v*,  v")  were  calculated  using  the  q(v*,  v“) 
values  with  r-ccntroid  correction.  The  solar  intensities  were  obtained  from  the 
Handbook  of  Geophysics  (I960). 

Curve  2.  Transition  strengths  as  in  Curve  I,  but  with  solar  intensities  equiva¬ 
lent  to  those  of  a  blackbody  at  6000*K. 

Curve  3.  Without  any  r -centroid  correction,  F  ■  0.  but  with  the  same  solar 
intensities  as  in  Curve  1. 

Curves  4  and  5.  With  the  intensity  ratios  computed  by  Authier  (1964)  and  Arm¬ 
strong  (196  3)  respectively,  where  the  transition  strengths  S(V,  v*)  were  calculated 
using  only  the  Franck -Condon  factors. 

From  the  five  curves  in  Figure  9,  the  following  conclusions  can  he  reached: 

1.  The  intensity  ratios  are  not  very  sensitive  to  errors  in  the  solar  intensities. 

2.  The  r -centroid  correction  is  important.  Without  this  correction,  the  de¬ 
duced  temperatures  will  be  far  too  high  (about  300 'K  for  the  -1  sequence  and  about 
100'K  for  the  0  Sequence). 


3.  There  are  some  discrepancies  between  our  computed  ratios  and  those  taken 
from  Authier  (1964)  and  Armstrong  (196  3).  The  latter  are  shown  as  dashed  curves. 
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9.  THE  ROTATIONAL  INTENSITY  DISTRIBUTION 


The  rate  of  absorption  into  an  upper  level  K'  is 


dNK' 

dt 


NxNy"  Nk"  ^  v(v,'K'>  v"«  K">  K't  V.K") 

=  NvMl)  Nv,i  S(v',  v")  ["n^i  ,  i  -Kill  v(v',K\  V",  K’ r  1)  0p  + 
x  3hcV'  [  K  +  l  2K'+  3  r 

+  NK'  -1  "<*,  K'<  v".  K'-l)eRj 

s  NxiSrA(v'  k,)’  where  nk"  =  -a^1  exP  •  W 


The  rate  of  emission  from  an  upper  level  K'  is 


dN 


KL 


dt 


=  NANV,  Nk,  K',v",  K")3  3(V,  V) 

K  ’  V" 

‘A  ^  ^,S(V''  '”>[«'  +  »  *'•  K’  f  >»3  + 

+  k'  v  (v1 ,  K\  v",  K'  -1)31  —  Na-^4nv,  Nk.  — E(v',  K’). 

J  A  3hc3  V  K  2K'  +  1 


=  N, 


(2) 


These  two  rates  are  equal,  and  the  population  in  the  upper  states  is,  therefore, 
given  by 


Na  N  ,,  Nk-  =  Nv  tel  (2K1  +  1)  , 

A  v  K  3hc3  X  3hc  E(v',  K’) 


(3) 


The  theoretical  intensities  are 


I(v',  K',  v",  K")  =  Nx  tel  S(V,  v">  S(K't  K">  K’-  v"’  K">3  (4) 

A  3hc  E  (v' ,  K  ) 


with 


iP  =  Nx  4f-3  Mvl*  £!• s(v’* v">  • (K' +  1)  •  v(v’>  K'* v"*  K’ _1)3 

v  A  3hc  E  (v',  K  ) 

■r  =Nx1nl  EMI' s' <v'-  -K'  •  *'V.  K’.  V.  K*  +  U3 


(5) 


(6) 
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In  order  to  compute  Ip  and  1^ ,  we  must  know  the  solar  intensities  corrected  for 
Fraunhofer  lines  for  transitions  which  are  strong  in  absorption. 

We  select  only  the  following  vibrational  transitions  which  have  to  be  corrected 
for  Fraunhofer  absorption:  v' -v"  =0-0,  0-1,  1-0,  and  1-1.  Thus,  about  500  lines 
of  importance  are  obtained. 


10.  BAND  PROFILES 

Using  Eq.  (4)  for  the  theoretical  intensities  of  the  lines,  the  band  profFej  for 
different  temperatures  can  be  computed. 

To  be  able  to  compare  these  band  profiles  with  the  recorded  ones,  a  convolu¬ 
tion  with  a  known  slit  function  is  performed.  The  empirical  slit  function  is  obtained 
by  letting  the  spectrometer  scan  a  single  line. 

Figure  10  shows  a  group  of  theoretical  spectra  for  temperatures  ranging  from 
300°K  to  900°K,  all  smeared  with  a  slit  function  of  2.  5  A  width.  The  spectral  range 
is  in  the  domain  of  Av  =  0,  with  three  vibrational  transitions. 

Figure  11  shows  band  profiles  from  the  Av  =  +1  sequence,  smeared  with  a  slit 

o 

function  of  4.  5  A  width.  The  band  profiles  were  corrected  for  Fraunhofer  absorption. 


Figure  10.  Synthetic  Spectrum  of  the  Av  =  0  Sequence  for  Various  Temperatures, 
2.  5  A  Slit  Function 
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Figure  11.  Synthetic  Spectrum  of  the  Av  =  +  1  Sequence  for  Various  Temperatures, 
4.  5  A  Slit  Function 


The  effect  of  Fraunhofer  absorption  on  the  emission  bands  is  demonstrated  in 
Figure  12  for  a  temperature  of  700°K.  These  profiles  were  obtained  by  using  the 
same  slit  function  as  for  Figure  10.  .  It  is  seen  that  the  two  profiles  differ  appreciably. 


II.  RADIATIVE  DYNAMICAL  EQUILIBRIUM 


The  deduction  of  temperatures  from  spectra  is  based  upon  the  assumption  that 
the  molecules  undergo  a  sufficient  number  of  collisions  to  thermalize  their  distribu-  ♦ 
tion  of  vibrational  levels  between  each  excitation -emission  process. 

To  see  what  the  vibrational  spectrum  would  look  like  were  this  not  the  case,  we 
assume  now  that  the  rate  of  collisional  redistribution  is  very  slow  compared  to  the 
rate  of  excitation  followed  by  emission. 

The  rates  of  the  electronic  transition  into  and  from  an  upper  level  are, 
respectively: 


8  TT3  y 

3hc 


and 


fdNA* 

I  dt 

*■  Jout 


jl  y! 

3hc3 


out 


Figure  12.  Effect  of  Fraunhofer  Correction  on  Synthetic  Spectrum 


These  two  rates  must  be  equal  for  stationary-state  conditions: 


n;  =  n.  §  . 

A  A  8ir?2  v 


Taking  into  account  the  vibrational  levels,  the  rate  of  transitions  into  an  upper 
level  v'  is  equal  to  the  rate  of  transitons  from  that  level: 


l 


N"  v.  S..  0.. 
J  IJ  iJ  tj 


=  Nl  it  Y  s..  V.  - 

i  v2  4*  ik  k 


Here 


vij  =  v(vi'  vj)»  s  vj)» 


and 


S..  =  S(v'  v?)  a  R  “  (1-0. 4<>  r(v'  v"))2  q  (v 


v"). 
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For  the  lower  level,  we  obtain  the  equation 


N'i  £  -ml  Sml  «ml  =  Nn  Snl  4  • 

By  eliminating  N*  and  N"  between  the  two  systems,  two  sets  of  homogeneous 
equations  are  obtained,  the  solution  of  which  gives  the  relative  vibrational  popula¬ 
tions  N[  and  Nj  in  the  upper  and  lower  electronic  state. 


where  W..  =  S..  v..  3  and  U  .  =  S  ,  v  .  0  . 

lk  lk  lk  ml  ml  ml  ml 


Let  us  consider  values  of  v*  and  v"  from  0  to  7.  The  two  systems  are  homo¬ 
geneous  sets  of  eight  equations  each,  and  it  enables  us  to  deduce  the  relative  values 
of  Nv,/Nv,  _0  a.id  Nv„  /NV„=Q  .  Figure  13  shows  the  values  of  Nv,/Nvl  sQ  plotted 
as  a  function  of  v*  (the  trace  marked  "8  levels"). 

Our  system  is  limited  to  eight  equations,  and  it  is  an  open  question  if  it  is 
sufficient.  In  order  to  see  the  trend,  we  reduce  the  system  to  six  and  then  to  four 
equations,  which  we  solve. 

Figure  13  shows  that  the  curves  marked  "6  levels"  and  "4  levels"  do  not  devi¬ 
ate  much  from  the  "8  level"  curve,  which  indicates  that  the  system  is  sufficiently 


large. 

Comparing  Figure  13  with  Figure  6,  it  is  seen  that  the  assumption  of  solar 
radiation  equilibrium  leads  to  relative  populations  in  the  upper  vibrational  levels 
which  are  far  greater  than  those  obtained  by  assuming  thermal  equilibrium  with 
the  ambient.  Since  the  number  of  collisions  required  for  vibrational  thermaliza- 
tion  is  usually  large  (10^  -  104),  it  would  appear  that  the  vibrational  temperature 
will  be  equal  to  ambient  temperature  at  altitudes  up  to  120  km,  where  the  collision 
frequency  is  about  150  see-1.  However,  experiments  performed  by  Blamont  et  ai 
(1961)  and  by  AFCRL  (to  be  published)  indicate  that  plausible  temperatures  can  be 
deduced  from  clouds  deposited  as  high  as  170  km  (approximately  one  collision  sec'*). 
This  means  that  most  of  the  upper  state  molecules  are  excited  from  thermalized 
ground  state  molecules.  This  could  be  due  to  continuous  formation  of  AlO  from 
the  released  TMA.  Clearly,  the  rate  processes  have  to  be  examined  in  more  detail 
for  accurate  temperature  measurements. 


i 
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Figure  13.  Relative  Population  for  the  Case  Where 
Lower  Levels  arc  not  Thermalized 


It  is  also  questionable  if  the  population  is  sensitive  to  errors  in  the  value  of 
the  transition  strength  S(v",  v").  By  leaving  out  the  r-centroid  correction,  the 
values  of  S(v',  v")  become  equal  to  the  Franck-Condon  factors.  The  solution  of 
the  set  of  equations  based  on  this  simplification  is  indicated  by  a  dashed  line  in 
Figure  13.  It  does  not  deviate  much  from  the  other  one. 

The  solar  intensity  changes  within  the  range  of  the  spectrum  of  AlO  by  a  factor 
of  ten.  To  determine  how  the  population  is  influenced  by  the  shape  of  the  solar  in¬ 
tensity  curve,  we  let  the  intensity  be  constant  and  obtained  the  curve  marked 
”0  (V,  v")  =  constant."  A  small  inaccuracy  in  the  solar  intensity  curve  would  not 
greatly  change  the  N'v,/Nvi  =  0  values. 
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